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ABSTRACT 

Flat radio spectra with large brightness temperatures at the core of AGN and X-ray binaries 
are usually interpreted as the partially self-absorbed bases of jet flows emitting synchrotron 
radiation. Here we extend previous models of jets propagating at large angles to our line of 
sight to self-consistently include the effects of energy losses of the relativistic electrons due 
to the synchrotron process itself and the adiabatic expansion of the jet flow. We also take into 
account energy gains through self-absorption. Two model classes are presented. The ballistic 
jet flows, with the jet material travelling along straight trajectories, and adiabatic jets. Despite 
the energy losses, both scenarios can result in fiat emission spectra, however, the adiabatic 
jets require a specific geometry. No re-acceleration process along the jet is needed for the 
electrons. We apply the models to observational data of the X-ray binary Cygnus X-l. Both 
models can be made consistent with the observations. The resulting ballistic jet is extremely 
narrow with a jet opening angle of only 5". Its energy transport rate is small compared to 
the time-averaged jet power and therefore suggests the presence of non-radiating protons in 
the jet flow. The adiabatic jets require a strong departure from energy equipartition between 
the magnetic field and the relativistic electrons. These models also imply a jet power two 
orders of magnitude higher than the Eddington limiting luminosity of a 10 M Q black hole. 
The models put strong constraints on the physical conditions in the jet flows on scales well 
below achievable resolution limits. 
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1 INTRODUCTION 

The centres or cores of many AGN show a flat radio spectrum in 
the sense that for the flux density as a function of frequency v we 
observe F v oc v a with a ~ 0. The high surface brightness tempera- 
ture associated with these spectra suggests a synchrotron origin of 
the emission. Observations with high spatial resolution reveals that 
the flat spectrum arises in the base of jet flows which continue to 
much larger scales (for a review see Cawtho rnelll99ll) . Similar flat 
or inverted (a > 0) radio spectra are also ob served in X-ray bina- 
ries in the low-hard state (e.g. Fendei, 2001). If optically thin, the 
flat synchrotron spectrum would imply a power-law energy distri- 
bution of the radiating relativistic electrons with a slope of unity. 
Such a distribution is very unlikely to arise for the usually assumed 
mechanism for the acceleration of the electrons at shock fronts (e.g. 
lBel ll978t) . 

A magnetised plasma containing very energetic electrons with 
a power-law energy distribution will produce a power-law spectrum 
at high frequencies. The slope of the spectrum, typically a < 0, is 
determined by the slope of the energy distribution. However, below 
a critical frequency the radiating electrons will re-absorb some of 
the photons. In this self-absorbed, optically thick regime the spec- 
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trum has a power-law slope of 5 /2, independent of the slope of the 
electron energy distribution (e.g. Rvbicki & Lightman. ll979l) . The 
spectrum of a uniform, self-absorbed synchrotron source therefore 
shows a pronounced peak. Blandford & Konial ( 1979) pointed out 
that in a jet the plasma conditions are changing along the flow and 
therefore the peaks of the self-absorbed spectra of different parts of 
the jet can occur at different frequencies. If the plasma conditions 
change such that the spectra peak at the same level, then the overall 
spectrum, observed with a spatial resolution insufficient to resolve 
the individual parts of the jet, will be flat. Their model has become 
the standard tool for interpreting observations of flat radio spectra 
from jetted sources. 

In the 1 Blandford & Konigll < 19791) model the jet is assumed 
to have a conical geometry, i.e. the velocity with which the jet is 
expanding sideways, is constant. The bulk velocity of the jet ma- 
terial along the jet axis is also assumed to be constant. The mag- 
netic field is assumed to be directed perpendicular to the jet axis 
and 'frozen' into the jet plasma. Adiabatic losses of the electrons 
are mentioned by the authors, but are assumed to be replenished 
by an unknown, continuous re-acceleration process along the en- 
tire jet. The same assumption is made for radiative energy losses 
associated with the emiss ion of synchrotron radiation. The subse- 
quent model of Marscher 1 1980) includes a simplified treatment of 
energy losses of the electrons due to adiabatic expansion and ra- 
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diative processes. It also allows for more confined jets, i.e. the jets 
are not necessarily conical. With these assumptions, the model is 
unable to pro duce a fla t emission spectrum. A similar model was 
developed by Hiellming & Johnston 1 1988). They consider adia- 
batic, but not radiative, energy losses of the electrons. The jet ge- 
ometry is again conical, but they also investigat e a more con fined 
jet. Th e model can predict flat spectra, but Hie llming & Johnstonl 
( 1988) point out that these may only arise under special circum- 
stances, particularly in the case of confined jets. The model of 
iGeorganopoulos & Marscher 1 1998) includes a detailed treatment 
of the energy losses of the relativistic electrons, but it concentrates 
only on the optically thin part of the spectrum of jets propagat- 
ing close to the line of sight for which numerical solutions are 
presented. The perhap s most com prehensive study of jet emission 
models is that o flRevnoldsliT98l) which includes the effects of en- 
ergy losses on the electron population, but neglects the effects of 
self-absorption on the electron energy spectrum. 

In this paper we extend the previous models by including adi- 
abatic and radiative energy losses and gains (due to absorption) for 
the electrons as well as investigating various possibilities for the 
evolution of the magnetic field. We consider two distinct cases: The 
ballistic and the adiabatic jet models. In the ballistic case the jet ma- 
terial follows straight trajectories and does not behave like a fluid, 
because individual fluid elements do not inte ract with each other. 
In ma ny ways this model is similar to the Blandford & Konigl 
1 19791) model, but we show that because of self-absorption effects 
we do not need to invoke a re-acceleration process to achieve flat 
emission spectra, in the adiabatic jet model the relativistic elec- 
trons suffer from adiabatic energy losses as well as radiative losses. 
Again we show that the models can produce flat spectra without 
re-acceleration of the electrons, but only for a very specific jet ge- 
ometry. The emphasis of our treatment is on the construction of an- 
alytical models and so we concentrate on jets propagating at large 
angles to the line of sight, i.e. the viewing angle is larger than the 
inverse of the Lorentz factor of the jet flow. 

In Section |2~T1 we briefly discuss the basic properties of our 
jets in terms of their geometry, the evolution of the magnetic field 
and that of the relativistic electrons. We present the first fully an- 
alytical solution of the equations governing partially self-absorbed 
synchrotron emission from a jet in Section l2~2l Section [23] sum- 
marises the model results for the case without radiative energy 
losses as studied in many previous models. In Section ISTTI we de- 
velop the formalism for including radiative energy losses in the 
model and the resulting spectra are discussed in Section 1X21 We 
apply the model to the data obtained for Cygnus X-l in Section 
[4] and derive the properties of this jet. Finally, we summarise our 
conclusions in Section|5| 



2.1 The basic jet properties 

2.1.1 Jet geometry and velocity 

We take the x-axis as the centre of a jet that is rotationally symmet- 
ric about this axis. The geometrical shape of the jet is then given 
by a one-dimensional function r(x) defining the jet radius with re- 
spect to the x-axis. Analogous to previous work we parameterize 
this function as r(x) = r$ (x/xq)"' = r^l" 1 , where xo is an arbitrary 
position along the x-axis defining the dimensionless coordinate I 
and ro is a constant scaling factor. The value of the exponent a\ 
depends on the details of the confinement of the jet. Confinement 
by ex ternal pressure is the simplest mechanism l Blan dford & Reesl 
1974) , but can lead to problems with the collimation of the jet 
iBege lman et allll984l) . Confinement by magnetic fields has also 
been suggested by various authors, but it is unlikely that magnetic 
fields alone, wi thout ad ditional gas pressure, can collimate the jet 
on large scales iBegelman, 1995). For our purposes here we do not 
need to specify the details of the jet confinement and we will as- 
sume that < a\ < 1. In principle one could also envisage highly 
overpressured jets with an accelerating expansion rate, i.e. a\ > 1. 
However, the pressure in such jets would fall very rapidly and they 
would quickly evolve to a situation where a\ < 1. 

The extreme case of a highly overpressured jet is that of a 
jet flow expanding into a (near) vacuum. In such a ballistic jet, as 
opposed to the adiabatic, confined jet discussed above, the jet freely 
expands in the direction perpendicular to the jet axis. In this process 
random, 'thermal' energy of the jet material is converted to ordered, 
kinetic energy associated with the sideways expansion. The random 
energy of the electrons giving rise to the synchrotron emission is 
reduced by this adiabatic expansion. However, in Section |4~T1 we 
apply the ballistic jet model to the observations of the jet in Cygnus 
X-l. There we will find a very small opening angle for the ballistic 
jet implying very small adiabatic expansion losses. Therefore we 
can assume that in the limiting case of a ballistic jet studied here 
the relativistic electrons do not suffer energy losses other than those 
associated with radiation processes. 

We assume in this paper that the velocity of the jet material 
along the jet axis, vj, is constant. While this is justified in the case 
of the ballistic jet, the confined, adiabatic jets can be accelerated, 
for example, by a pre ssure gradient in the confining medium. In 

the model of Blandfo rd & Reesl Il974t) the Lorentz factor of the 

-1/4 

bulk velocity is proportional top* , where p x is the pressure of 
the external medium. As long as the external pressure gradient is 
shallow, the Lorentz factor of the jet flow will be only a very weak 
function of the position along the jet axis. Similar arguments hold 
for a magnetically confined jet. The constant bulk velocity of the 
jet also implies that a given volume element AV travelling with 
the jet flow will only expand sideways according to AV °= r 2 . A 
constant jet velocity also simplifies the model greatly as we can 
ignore the effects of varying length contraction along the jet axis 
iGeorganopoulos & Kazanaa.120041) . 



2 THE MODEL 

In this Section we derive the emission properties of partially self- 
absorbed jets neglecting radiative energy losses of the relativistic 
electrons. Note that we are concentrating on jets at comparatively 
large angles to our line of sight, f). As we will point out further 
down this greatly simplifies the determination of the optical depth 
of the jet material. 



2.1.2 Magnetic-field 

The strength of the magnetic field changes during the sideways ex- 
pansion of the jet material. In general, we parameterize the evolu- 
tion of the magnetic field as B = BqI~" 2 . For flux freezing of the 
magnetic field and using flux conservation we have that the field 
component parallel to the jet axis, Bn , is proportional to r~ 2 . Also, 
the magnetic field components perpendicular to the jet axis, B ± , 
are proportional to r~ 1 . For an initially mixed field, B^ will always 
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become the dominant component and so 02 = a\ . The perpendicu- 
lar magnetic field may also contribute to the confinement of the jet. 
For completeness we also consider a purely parallel configuration 
of the magnetic field with 02 = 2a 1 . Finally, if the magnetic field is 
constantly tangled by turbulent motions in the jet material on scales 
smaller than the jet radius, then it can remain isotropic and it be- 
have s like a relativistic fluid w ith B = BqI^ 4 " 1 ^ and 02 = 4ai/3 
(e.g. lHeinz"& Begelman, 2000). This is analogous to the behaviou r 
of the magnetic field in an isotropic expansion (e.g. Longair, 1994), 
but is clearly incompatible with flux freezing. The last case of a per- 
manently isotropic field cannot be realised in the ballistic jet as it 
would require that the jet material behaves like a fluid. 
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Table 1. Exponents of the model parameters used in this paper. See text for 
details. 



2.1.3 Relativistic electrons 

In order to produce synchrotron emission the jets must contain a 
population of relativistic electrons. We assume that the latter has a 
power-law energy distribution of the form 



N(E)dE = KE~ p dE, 



(1) 



where E is the electron energy, E = "fln^cr, and K is a scaling in- 
dependent of E. y is the Lorentz factor associated with the rela- 
tivistic motion of the electrons. In this section we do not impose a 
high-energy cut-off on the energy distribution and we neglect radia- 
tive energy losses. Even so the energy distribution of the electrons 
changes as the jet expands. We represent the evolution of the elec- 
tron distribution by setting K = KqI~" 3 . 

For a given volume of jet material AV particle conservation 
demands that 



AVKy p dy = AV KoY(/ d Yo> 



(2) 



where all quantities with subscript '0' refer to their values at x = xq . 
Therefore for the ballistic jet we have aj, = 2a\ =2. 

For the adiabatic jet we need to include energy losses due to 
the jet expan sion. Since mos t of the electrons are highly relativistic 
we have (e.g. lLongaiilll994t) 



dj 1 31b (AV) 



dt 3 T 



dt 



which has the solution 

-1/3 



Y = Yo 



/AV_Y 



and it follows that 

37b = /AVV /3 
3y \AVq) ' 

Re-arranging equation J2j and substituting yields 



_. . AV fAV\- p/i _.( 



AV 



1/3 



dy. 



(3) 



(4) 



(5) 



(6) 



Collecting terms and remembering that AV 0= r{x) 2 we find 03 = 
(4+2p)oi/3. 



2.1.4 Individual models 

On the basis of the discussion above we formulate five individual 
models distinguished by the magnetic field behaviour. The ballistic 
jet models, B 1 and B2, as well as the adiabatic models, Al and A2, 
correspond to a perpendicular and parallel field structure, respec- 
tively. The adiabatic model A3 represents the case of an isotropic 



magnetic field in the jet. The relevant coefficients describing the jet 
geometry and the behaviour of the magnetic field and relativistic 
particles are summarised in TableQ 



2.2 Partially self-absorbed synchrotron emission from a jet 

From the expressions for r, B and K defined in the previous Section, 
we can now build a model for the emission from the jet. For this 
purpose we split the jet into small segments of length dx along the 
jr-axis. We assume that the segments move along the jet axis at a 
constant velocity vj = fijc corresponding to a Lorentz factor Yj- The 
jet axis is at an angle ■& to the line of sight of the observer and 
so the Doppler factor for an approaching (' — ') or receding ('+') 
jet is 8 T = [Yj (l TPjCos 1 ?)] ■ The observable monochromatic 
intensity of one such segment taking into account absorption is 



■/v 



47tXv 



1-e 



-Zvto 



(7) 



where J v is the emissivity per unit volume and %v is the absorption 
coefficient. For convenience in the development of the model the 
frequency v is measured in the restframe of the jet material. It is 
related to the observing frequency by v Q t, = 8 T V. Here and in the 
following we assume that the average path of a photon through the 
jet has the length r(x). An exact calculation of the radiative transfer 
of photons through various jet elements would have to take into 
account relativistic aberration effects. It is therefore complex and 
impossible in an analytical model. The assumption of an average 
path length r will not introduce a large error as long as the angle to 
the observer's line of sight is large. 

The jet segment has a surface area of 2nr(x) dx and so the 
observable flux density of the segment is given by 



dFv = 5 3ZWlYfl_ e -ZvrW 
T 2£> 2 Xv 



dx. 



(8) 



where D is the distance of the jet from the observer. 

Substituting the dimensi onless variable / = x/xo, we can ex- 
press J v and Xv in SI units as iLongaut [l994T) 



j v _ j v {l-p)/2j-a 3 -a 2 (p+l)/1 
^ v _ j^ )V (-p-4)/2f-a<,-o 2 (p+2)/2 



(9) 
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^0 

Xo 
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(10) 
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and the constants c\(p) and ci(p) given by equations 18.49 and 
18.74 in lLoneaiil fi994). Substituting into equation {SJ and integra- 
tion gives the total flux density of the jet as 



S 2 x Q r 0^0 



v 5/2 ^ iai+ a 2 /2^_ e -^ 



d/, 



(11) 



where / m ; n and / max are the physical limits of the jet flow along the 
x-axis and the optical depth of the jet material is given by 



t(0 = JCvOO = Xcov 



{-p-4) /2 [a , -a 3 -a 2 (p+2) /2 _ 



(12) 



The reduction in the number of Doppler factors arises from our 
assumption of a steady state of the jet flow. In principle further rel- 
ativistic corrections must be applied in the case of mixed optically 
thin a nd thick emission, but these corrections are small I Cawthorne, 
Il99ll) and we neglect them here for simplicity. 

It is convenient to recast equation 1 1 1 > with the help of equa- 
tion 1121 as an integration over optical depth, 



where 



A max 



Ms — 1 



1-e" 



dx, 



a\ = 1a\ — lay — (p + 2) 02 
and 
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(13) 



(14) 



(15) 



The coefficients 04 and ag are listed for the ballistic and adiabatic 
jets in Table Q to is given by setting / = 1 in equation 1 12i while 
Xmax = x(/ mm ) ar, d ^min = x(Z m a X ), which reflects the fact that the 
optical depth is always greatest in the innermost regions of the jet. 
The solution of equation 1 1 3i is given by 
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Here, the incomplete T-function is defined as 
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(17) 



2.3 Spectra in the absence of radiative energy losses and 
without a high-energy cut-off 

We can immediately recover the well-known solutions for an en- 
tirely optically thin (x mm <C 1 and x max C 1) and an entirely op- 
tically thick (t m i n 2> 1 and x max ^> 1) jet. We note that for all 
choices of a\ discussed above and for physical reasonable val- 
ues for the power-law exponent 2 < p < 3 we find ag < (see 
Table 0. The incomplete T-function has the series representation 
I Gradshtevn & Rvzhikl EoOO) 



(a.z) 



-r(a) 



~ n\{a+n) 
For x-C 1 we ignore all terms beyond n ■ 



1 and thus obtain 
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(18) 



(19) 
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Because of equation J 121 we have Xo °c X^,, X max 0= 
and from equation I16i it then follows that F v oc yi^-p)/ 2 as ex- 
pected. 

For large optical depths we ca n use th e limit for the incom- 
plete T-function, lim^^, r (a, z) = 1 Grads htevn & Rvzhild. ,2000). 
Thus 
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Figure 1. Relation between the geometrical shape of the jet of model A3 
and the power-law exponent of the energy distribution of the relativistic 
electrons for which the jet emission spectrum is flat. The plot shows the 
relation in equation 1241 . 
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and from the proportionality of the optical depths terms it then fol- 
lows that F v <x v 5//2 , again as expected. 

The final special case is that of a spatially very extended jet or 
'long' jet. If the physical dimensions of the long jet, / m j n and / max , 
are such that x max — * 00 and x mm — » 0, then for —1 < ag < we 
have 



r( fl5 ,x) + -x fl5 

ag 

which implies 

F v oc V [ 5 +("+ 4 ^]/ 2 . 



' r («5), 



(21) 



(22) 



Using the results summaris ed in Table Q we recover the result of 
Blandford & Konigil 1 1979b that the ballistic jet with the magnetic 
field perpendicular to the jet axis (model Bl) has a flat spectrum, 
i.e. F v is independent of v, if it is extended and — 1 < ag < 0. The 
ballistic jet with a parallel magnetic field, model B2, can never pro- 
duce a flat spectrum for positive p because equation J22> predicts 
F v oc \0p+' } )/1 2 ( 2 p+' } )}. For the adiabatic jet models we can substi- 
tute the expressions for ag and find that a flat spectrum is predicted 
if the exponent a for the geometrical shape of the jet, r(x) « x°, is 
given by 
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(23) 



15p + 5 

FigureQplots the relation for model A3. It is interesting that for all 
adiabatic jets geometrical shapes described by exponents a in the 
range l/3;$a<2/3 are required for flat spectra. 



Note that for ag < 
J 1 8i dominates for x m ; n 
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1 the n = 1 term in the series in equation 
-> 0. In that case we have for the long jet 
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Figure 2. Flux density per unit length (solid line) at a single frequency as 
a function of position along the jet axis for the ballistic jet. This plots the 
integrand in equation 1111 . The optical depth at the same frequency is also 
given (dashed line). 
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Figure 3. Spectrum of the ballistic jet model B 1 . The solid line shows the 
spectrum without radiative losses or high-energy cut-off in the electron en- 
ergy density distribution. The short-dashed line is the same spectrum in- 
cluding synchrotron losses and a high-energy cut-off. The dotted lines show 
the frequency at which the optical depth at the beginning (high frequency) 
and the end (low frequency) of the jet, / ra j n and / raax respectively, is equal 
to unity. The dot-dashed line at low frequencies indicates where the optical 
depth is equal to 10~ 4 at / max . The dot-dashed line at high frequency shows 
where the optical depth is 100 at / m ;„. 

similar to the entirely optically thin jet. The spectrum of the long 
jet is then also optically thin, i.e. F v oc y^^P)/ 2 , and a fiat spectrum 
is not possible for physically reasonable values of the exponent p. 

Unless the jet is exceedingly short, there will always be a 
range of frequencies for which x m j n — > and t max — > oo and the 
long jet scenario applies. An example is shown in Figure|2|where 
we plot the integrand in equation < I I > and the optical depth of the 
jet material as a function of / for a single frequency for the ballis- 
tic jet. The contribution to the overall flux of the jet peaks close to 
T = 1. For the long jet scenario to apply the jet must be long enough 
so that substantial emission from either side of the peak contributes 
to the overall flux. 

Below the frequency range of the long jet the spectrum will 
follow the optically thick case, F v v 5//2 , and above this range the 
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Model parameter 


Value 


xo 


47 AU 


ro 


8.9 x 10 7 m 


P 


2.5 


B 


2.4 mT 


KO 


3.3 x 10- 7 J 15 m 


D 


2 kpc 


An in 


4.3 x ur 5 




200 


Ymax (^min) 


10 6 


v i 


0.97 c 


6 


40° 



Table 2. Parameters used to illustrate the spectral properties of the jet mod- 
els. These model parameters are also used to explain the observational data 
of Cygnus X-l in Section l4*Tl when using a ballistic jet model with a mag- 
netic field perpendicular to the jet axis, model B 1 . 



1 000.0 F 




Frequency / Hz 

Figure 4. Same as Figure |3] but for the adiabatic jet model A3 with a = 
0.46. The long-dashed line shows the spectrum including both adiabatic 
and synchrotron energy losses. 



optically thin case applies with F v ac vi^i^/ 2 . The solid line in 
Figure [3] illustrates this generic overall shape of the jet spectrum 
for the ballistic jet model Bl. For this Figure and the following 
we have used the model parameters summarised in Table [2] The 
observations of the partially self-absorbed jet of Cyg X-l are well 
explained by the model for these parameters (see Sections l3.2l and 
|4~TY For comparison, the solid line in Figure|4]shows the spectrum 
of the adiabatic jet model A3 with a = 0.46 for the same set of 
parameters. The value of a was chosen according to equation 122\ 
to allow for a flat spectrum at intermediate frequencies. 

It is interesting to note that the frequencies for which the op- 
tical depth of the jet material is unity at l m [ n and / max are located 
well within the optically thick and thin regimes respectively. The 
transition from the, in this case, flat spectrum of the long jet occurs 
closer to t(Z max ) ~ 10~ 4 and x(/ m i n ) ~ 100. Obviously, for suffi- 
ciently short jets the frequency range over which the long jet case 
applies may vanish altogether. 
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3 INCLUDING ENERGY LOSSES OF THE ELECTRONS 

In the previous Section we did not consider the effect of radiative 
energy losses of the relativistic electrons on the predicted spec- 
tra. The adiabatic jet models include the effect of adiabatic energy 
losses on the overall energy distribution of the relativistic electrons. 
However, because we did not impose a high-energy cut-off to this 
distribution, we did not have to consider the effect of adiabatic 
losses on such a cut-off. In this Section we introduce a high-energy 
cut-off at Ymax and include the effect of adiabatic energy losses on 
this cut-off. 



3.1 Evolution of the high-energy cut-off 

3.1.1 Adiabatic and synchrotron losses 

Other than adiabatic energy losses, the radiative losses due to syn- 
chrotron radiation modify the energy distribution of the relativis- 
tic electrons awa y from a simple power-law, unless p = 2 (e.g. 
iKardashevl Il962l) . In the following we will make the simplify- 
ing assumption that the energy losses only shift the sharp high- 
energy cut-off to lower energies while not altering the power-law 
shape or exponent of the power-law distribution. This approxima- 
tion does not introduce a large error as the deviation from the orig- 
inal power-law is significant only near the cut-off. Also, the ex- 
pressions for the synchrotron emissivity and absorption coefficient 
given in equations {5J and J 101 are strictly valid only for power-law 
energy distribution extending from y = 1 to Ymax — » °°. However, 
the expressions involved in the derivation of J v and % v decay suf- 
ficiently quickly for ysi 10 that the results for finite Ymax do not 
deviate gr eatly from those presented in the previous Section for 
Ymax -> °° iRvbicki & Lightmarl[l97^lLongaiil.ll994l) . 

In the optically thin regime the evolution of the Lorentz factor 
of a given relativi stic electron in the rest frame of the jet material is 
described by (e.g. lLongaiil.ll994t) 



4 Ox i( 
3 m e c Wq 



-2a-i 



2a i 
~3f 



(25) 



where the first term on the right describes the energy losses due 
to synchrotron radiation and the second term reproduces equation 
for the adiabatic losses where we substituted for AV. Because of 
our assumption of a constant bulk velocity for the jet material along 
the x-axis, vj, we can express the dimensionless coordinate / = x/xq 
also as a time variable, i.e. / = YjVj'Ao- Because of time dilation, 
we have to include Yj- Thus t = Ixq/ (VjYjJ and?o=-W ( v jYj)-°Ti s 
the Thomson cross section and uq = Ek/ (2/jq) is the energy density 
of the magnetic field at xo. The solution of equation I25i is found 
.i 



y(0 



Y('min)f- 2fll/3 



-2o,/3 tori in , \(m-t a( ' )' 
r min + 3fl 6 m c e r H r minJ (» r min ) 



(26) 



with «6 = 1 — 2«2 — 2ai/3 and t mm = / m in-W ( v jYj)- Electrons 
which were injected into the jet at t mm or, equivalently, / m [ n with 
a Lorentz factor Y( f min)> have a Lorentz factor y(f) at t or, equiva- 
lently, /. For the ballistic jet the adiabatic, second term on the right 
of equation 1251 vanishes and we have instead 
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Note that the exponents ag and 1 — 2«2 are usually negative. 
This implies in the case of the ballistic jet models (no adiabatic 



losses) that the Lorentz factors of electrons do not necessarily de- 
crease forever, but converge to a finite value for t f m j n . The some- 
what surprising result simply reflects the fact that the synchrotron 
losses rapidly decline in the decreasing magnetic field of the ex- 
panding jet. For the adiabatic case the adiabatic losses continue at 
all times and so Y°= f~ 2fl '/ 3 oc /- 2fl i/ 3 for t > t^. 

For optically thin conditions the evolution of the high-energy 
cut-off Ymax also obeys equations 1261 and 1271 . Below we will re- 
fer to the high-energy cut-off in the optically thin regime as Ythin- 
However, for large parts of the spectrum the jet is optically thick. 
Electrons with a given Lorentz factor y emit most of their radia- 
tion at the critical frequency v ~ VgY 2 , where the gyro-frequency 
is defined as v g = eBj (2nm e ). An electron emitting at a critical 
frequency for which the jet is optically thick gains energy through 
synchrotron self-absorption. Ideally we would include an energy 
gain term for the self-absorption effect into equation J25t and then 
derive the electron evolution as before. While this approach leads 
to analytic solutions when only considering the systematic energy 
gain of electrons of a single energy lReest ll967). it is not applica- 
ble in most cases because the stochastic energy gain for power-law 
energy distributions is comparable to the systematic term. In this 
case, only numerical solutions are possible, be cause the stoch astic 
term depends on the entire energy distribution iMcCravll 19691) . 

The full numerical treatment of synchrotron losses and gains 
in the optically thick regime is beyond the scope of this paper. How- 
ever, electrons radiating mainly at frequencies for which the jet is 
optically thick, do on average not lose or gain energy due to radia- 
tive effects, eve n if they are relatively close to the surface of the jet 
jMcCravl ll969). Thus the high-energy cut-off in the optically thick 
regime, Ythick* ls given by the requirement that T(ythick) ~ L The 
electrons at this cut-off emit mainly at a frequency V tmc k = v gYthi c k 
and so we find from equation < 1 21 
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For the ballistic jet with perpendicular magnetic field, model Bl, 
a-j — and so Ythick ls constant along the entire length of the jet. 
In other words, electrons with Lorentz factors equal or below Ythick 
never loose their energy to radiation unless they are very close to 
the jet surface. The existence of a constant hig h-energy cut-off is re- 
quired for a flat spectrum from the jet. While Blandford & Konigl 
1 1979) invoked an unknown re-acceleration mechanism to ensure 
Ymax = constant, we have shown here that such a process is un- 
necessary because of the energy gains associated with synchrotron 
self-absorption. For a parallel magnetic field in the ballistic jet, 
model B2, we have a-j = 2/ (p + 4) . The Lorentz factor of relativis- 
tic electrons for which the jet is optically thick is increasing for 
increasing / in this model. Therefore, if / m i c k is the position along 
the jet where Ythin = Ythick, we have Ymax = Ythick at this position 
and y max = constant for all / > / t hkk- In me adiabatic cases we find 
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Even for optically thick conditions the electron energy distri- 
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bution does not deviate greatly from the original power-law with a 
high-energy cut-off y max for exponents 2 < p < 3 <McCravl[l969l) . 
Thus, for each position / along the jet we can now determine Ymax 
and thereby the entire electron energy distribution. For small Z the 
cut-off is given by Ymax = Ythin- Further down the jet Ythin will nrs t 
become equal to and then fall below Ythick an d f° r the ballistic jet 
models Ymax = Ythick afterwards. In the adiabatic jet models two 
competing effects can diminish Ymax further after passing through 
the point at which Ythick = Ythin- m most cases equation I28i im- 
plies a further reduction of Ythick f° r increasing /. This means that 
Ymax would also decrease. At the same time, adiabatic losses lead to 
Ymax = Ythick ('/'thick) ~ 2a '^ 3 - K is straightforward to show that for 
all adiabatic models the adiabatic losses of Ymax are the dominant 
effect. 

By again making the assumption that all electrons only emit 
at their critical frequency, we can now define for a given frequency 
V a maximum distance / max along the jet axis where the jet material 
is still contributing to the overall emission, 

f 2nm e v V 1/(12 

(31) 
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Clearly, as long as / max is larger than the physical extent of the 
jet, / m ax> the jet spectrum is not affected by energy losses of the 
relativistic electrons at frequency v and we can use the results of 
the previous Section. For / max < / max we have to take into account 
energy losses of the electrons by using / max instead of Z max in the 
calculation of z m [ n . 

3.1.2 Losses due to Compton scattering 

In the optical thick parts of the jet Compton scattering of the syn- 
chrotron photons off the relativistic electrons may become impor- 
tant. We do not include energy losses of the relativistic electrons 
due to Compton scattering in the jet in our calculations as these 
would require a full treatment of radiative transfer. However, it is 
obviously necessary to test whether these losses are important when 
applying the model to observational data and so we give the neces- 
sary expressions below. 

The relevant limit for the energy density of the synchrotron 
photon field is most conveniently expressed in terms of the bright- 
ness temperature (e.g. lRvbicki & Lightmanlfl979T) . 
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where k B is the Boltzman constant and v is the emitted frequency 
rather than the observing frequency. For 7j-,jSl0 12 K Compton 
losses are not important compared to the energy losses due to 
synchrotron radiation. The maximum brightness temperature for a 
given frequency is reached at the position along the jet where the 
optical depth of the jet material roughly equals unity for photons of 
this frequency. Hence in our model we have 
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The maximum brightness temperature is either constant or only a 
weak function of I in all our models. Also, using the dependencies 
of Jq and Xo given in equation 1 101 we find that 7}, max depends only 
weakly on the other model parameters Bq, Kq and tq. 

In some jets relativistic induced Compton scattering may be 
more important than direct Compton scattering discussed above. 
Induced Compton scattering causes signific ant en ergy losses for 
the relativistic electrons if iSincell & Kroliklll994) 
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where tj is the Thomson depth of the jet material, 
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Here, n e is the number density of electrons, o"t is the Thomson 
cross-section and we have again assumed that the average path 
length a photon travels through the jet material is equal to the jet 
radius, r. For our power-law energy distribution of the electrons an 
upper limit for the electron density is given by n e < K (m e c 2 ) . 
The Thomson depth of the jet materials in our models is then lim- 
ited by 

1-P 



O T K r r 



(36) 



Equations I33i and J36I can be used to ensure that the models 
are applicable to a given observational data set, i.e. that they do not 
suffer from Compton losses which are not included in the models. 

3.2 Spectra with energy losses and a high-energy cut-off 

3.2.1 Iterative construction of model spectra 

For a given set of model parameters we can construct a model spec- 
trum. In practice this will involve the determination of T mm and 
Tmax from equation J 1 21 to be substituted into equation 1161 . For a 
given frequency v, we set / = / m i n and calculate T m ax- Determin- 
ing T m j n is more involved as it requires the calculation of / max . This 
calculation involves an implicit equation and so cannot be done an- 
alytically. Here we describe one possible iterative procedure for the 
determination of / max . 

The first step is to choose a trial distance Z such that < 
' < 'max- The strength of the magnetic field in the jet material at 
I is B(l) = B()l~ ai . The maximum frequency at which jet mate- 
rial located at / is still contributing to the emission is given by 
Vmax(') = Vg (I) Ymax (')• F° r the next iteration we need to com- 
pare v max (/) with v. Therefore we must next derive Ymax (')> the 
high-energy cut-off of the electron energy distribution at I. 

From equations J26I adiabatic jet models) or J27I ballistic jet 
models) we can determine Ythin (')• We calculate Ythick from equa- 
tion J28j- If Ythin (0 > Ythick, then Ymax (') = Ythin (')• Otherwise, for 
the ballistic jet models Ymax (') = Ythick- F° r the adiabatic jet models 
Ymax (') = Ythick ('/ 'thick) -2 " 1 ' 3 • The required distance l± ick must be 
found from the implicit equation resulting from setting j(t) = Ythick 
in equation 1261 . Finally, if V m ax (') > V, then the trial distance in 
the next iteration should be larger than the current one. In the case 
of v max (/) < v, the trial distance should be decreased. The iter- 
ations can be stopped when v max (') ~ V within the required ac- 
curacy and at that point we can set l' m . dx = / and then proceed to 
calculate T mm . 



3.2.2 Example spectra 

The model parameters in Table [2] were chosen to explain the ob- 
servations of the jet in Cygnus X-l with the ballistic jet model Bl 
including radiative energy losses (see Section |4~T1 . This does not 
imply that energy losses will always be important in all jets and at 
all frequencies. 

The spectrum of the ballistic and adiabatic jets with energy 
losses of the electrons are shown in Figures[3]for model Bl and|4] 
for model A3. When energy losses of the electrons are taken into 
account, then we cannot in general expect that z m [ n <C 1 for a given 
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frequency. Therefore the spectrum of the jet will not necessarily be 
that of the long jet described by equation I22i . In fact, in most cases 
T m i n will considerably exceed unity. If T max — > °°, then equation 
\2()\ applies and we find 
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to 



(37) 



For the ballistic jet models we have argued in the previous 
Section that Ymax = constant for / > /thick- F° r rapid radiative energy 
losses of the electrons at / in the range / mm < / < /thick- the distance 
/thick itself will not depend strongly on the observing frequency. 
From equation J3 1 1 we then find / max ^ v -1 /" 2 and substituting 
into equation 1121 we get 

Vm-v^- 1 "^)/" 2 . (38) 

Note that the exponent does not depend on the slope of the elec- 
tron energy distribution, p, nor on the geometrical shape of the jet 
described by a (see Table0. Finally, from equation 1371 we obtain 
the slope of the spectrum as 

(4a 2 -2ai-2)/(2a 2 ) 



(39) 



For the ballistic jet with perpendicular magnetic field, model Bl, 
we have F v = constant as in the case without energy losses of the 
electrons, which is confirmed by Figure|5| For a parallel field struc- 
ture, model B2, the spectrum would follow F v <^ v. Model B2 is 
incompatible with a flat spectrum. The slope in the optically thin 
part of the spectrum in Figure|5|is steeper compared to the case of 
no energy losses because of the effect of the high-energy cut-off in 
the electron energy spectrum. 

In the case of the adiabatic jet models y max ^ 

r 2a,/3 an( j 

SO 

'max x v~ 3// ' 4ai+3fl2 '. Again substituting into equation 1121 yields 



. v -(7ai+3a 2 )/(4fli+3a 2 ) 



(40) 



where we also used «3 = (4 + 2p)ai/3 as appropriate for the adi- 
abatic jet models. Again the exponent of this expression does not 
depend on p or a\ . The shape of the spectrum is now 

F v oc v (7ai+6a 2 -3)/(4<!i+3<i 2 )^ 

The exponent of the power-law spectrum predicted by model A3 
for ci\ = 0.46 is then 1.06 which is confirmed by the slope of the 
spectrum in Figure [4] below about 10 12 Hz. The emission at high 
frequencies comes from the innermost parts of the jet close to Z m j n . 
Radiative losses had not enough time there to completely change 
the electron energy distribution. This explains the peak in the spec- 
trum. At the highest frequencies only optically thin parts of the jet 
contribute to the overall emission and lead to a negative power-law 
similar to the case without energy losses. Note however that the 
slope of this power-law is somewhat steeper due to the decreasing 
high-energy cut-off. 

The adiabatic jet models are all consistent with a flat spectrum 
provided the shape of the jet described by a\ takes a suitable value 
(Model Al: a x = 3/13, A2: a x = 3/19 and A3: a x = 1/5). In all 
three adiabatic models the jet needs to be strongly confined, i.e. 
a\ < 1/4, to achieve a flat emission spectrum. The slope of the 
spectrum is quite sensitive to the value of a\ . For example, a change 
of a\ in model A3 from 1 /5 to 1/3 results in a change of the power 
law exponent of the spectrum from zero to 0.75. 

For the following discussion we note that all the spectral 
slopes calculated above, with the only exception of the purely opti- 
cally thin case, are all independent of the exponent p of the power- 
law describing the energy distribution of the electrons. Therefore 
we can readily apply the model to observational data even if the 



optically thin part of the spectrum is not observed and thus we do 
not know the value of p. 



4 APPLICATION TO OBSERVATIONS 

The jet emission models depend on a number of parameters. Some 
of these parameters can be constrained by applying general con- 
siderations and others may be inferred from applying the model 
predictions to observational data with a view to determining the 
physical conditions within the jet. Here we discuss all of the rel- 
evant parameters in turn and demonstrate below how observations 
of the jet in Cygnus X-l may be used to infer the properties of this 
object. 

The scale height xq can always be chosen arbitrarily to pro- 
vide a convenient location along the jet axis at which to define the 
exact values of other quantities. In many cases we will be mainly 
interested in that part of the jet spectrum which is strongly affected 
by absorption. As we have seen in the previous Section, we then 
do not need to know the exponent of the power-law energy distri- 
bution of the electrons, p, as it does not influence the slope of the 
predicted spectrum. However, estimates for other quantities derived 
from the model, for example the strength of the magnetic field, de- 
pend weakly on p. The distance of the jet, D, the bulk velocity of 
the jet material, v,-, and the viewing angle of the jet axis to our line 
of sight, r>, cannot normally be determined by the model itself and 
need to be measured by other means. 

The parameters describing the geometrical shape and size of 
the jet, tq, Z m j n and / max , could in principle be determined from 
observations. However, / m ; n is probably too small to be resolvable 
even with a large improvement on current resolution limits. Cur- 
rently only upper limits exist for the radius of jets in Galactic X-ray 
binaries (Miller- Jones et al., in preparation) while for AGN jets rrj 
is sometimes resolved fe.g. ljunor et allfl999h . The maximum ex- 
tent of a jet at a given observing frequency is sometimes measured 
and an example i s provided by the observations of the jet in Cygnus 
X-l o flStirling et alJ l2001 ) which we use in the following Section. 
It should be borne in mind that at one observing frequency we can 
always only measure l' m3X given by equation 13 li rather than the 
physical extent of the jet flow / max . However, / max only determines 
the low frequency cut-off of the spectrum, but is not important for 
the model otherwise. In the case of the adiabatic jet models resolved 
observations can, in principle, also determine the shape of the jet 
as described by the parameter a. However, in practice it is easier 
to infer the value of a from the slope of the observed self-absorbed 
spectrum as this is a strong function of a. 

Finally, the normalization of the electron energy distribution, 
Ko, and the strength of the magnetic field, Bo, cannot be determined 
directly from observations, but must be inferred from the model. 
We can reduce the number of free parameters by assuming that the 
energy densities of the magnetic field and of the relativistic elec- 
trons are initially in equipartition. In this case, 
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For an energy distribution extending over all physically meaningful 
Lorentz factors (1 < y < °°) we then have 

2 

(43) 
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With these considerations we can now apply the model to ob- 
servations and determine relevant parameters for the observed jet. 
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4.1 Application to Cygnus X-l 

In the following we app ly the model to the jet observed in the X- 
ray binary Cygnus X-l. lStirling et alj 1200 lh report a resolved jet 
extending to about 15 mas from the position of the X-ray binary 
system at an observing frequency of 8.4 GHz. We set the bulk ve- 
locity of the jet material to 0.97 c and th e viewing angle to our line 
of sight to r> = 40° istirling et all EoOll) . For a distance of 2kpc 
( Gierl ihski et allll999l) the observed, projected jet length then cor- 
responds to a real jet length of roughly 47 AU. For convenience we 
set xq equal to this value and so / max = 1 at 8.4 GHz. Only one jet 
is observed and it is therefore reasonable to assume that this is the 
approaching jet. Note that the bulk velocity of the jet and the view- 
ing angle imply 8_ ~ 1 and sov~ v b- The total flux density at the 
same frequency is 13 mjy. 

In a map at 15 GHz from an earlier observing epoch, 
the jet may also be marg inally extended along the same axis 
( Stirl ing et all 1 1 998L EoOlt) . The extension is 2 mas or less which 
corresponds to / max (15GHz) < 0.13. In the following we will 
mainly concentrate on the observational data at 8.4 GHz. 

There are no simultaneous observations at any frequency other 
than 8.4 GHz and so we cannot be certain what the spectral slope 
of the jet emission was. However, during the low/hard X-ray state 
the source usual ly shows a fl at spectrum extending up to at least 
220 GHz jFender et all l200fj) . We assume here that the spectral 
slope at the time the radio jet at 8.4 GHz was observed was zero. 
The extent of the flat spectrum to high frequencies is also not 
known. However, for G X 339-4 the flat jet spect rum is observed 
to near-IR wavelengths l Corbel & Fender, 2002) while for XTE 
Jl 118+480 it may extend to the near-UV favnes et aill2000l) . For 
the purpose of illustrating the model, we assume here that the flat 
spectrum extends to the near-IR of wavelengths of about 1 pm. The 
flat spectrum always arises from regions in the jet which are at least 
partially self-absorbed. Thus the exact value of p is not very impor- 
tant and we set p = 2.5. 



4.1.1 Ballistic jet models 

We have seen above that the ballistic jet with a parallel magnetic 
field configuration, model B2, is inconsistent with a flat spectrum. 
In this Section we therefore concentrate on model Bl, a ballistic 
jet with a magnetic field perpendicular to the jet axis. With the as- 
sumptions made above equation I37i reduces to 

F v ~ 5.9 x 10- 7 5|x r Z)- 2 B 1/2 v 5 / 2 4 ax 5/2 mJy (44) 

for the ballistic jet model Bl. Here and in the following, all quanti- 
ties are measured in SI units unless indicated otherwise. Substitut- 
ing the measurements discussed above we get an expression for r$ 
as a function of By. A second equation relating the same quantities 
can be found from substituting Ythick from equation I28i for 
equation <3 1 1 resulting in 

C = 2.0xl0 11 V- 1 Br^ /13 '-0 /13 = l J (45) 

where we again made use of the observed quantities. For initial 
equipartition we can eliminate Kq and solve for Bq. We can then 
calculate all other model parameters and they are summarised in 
Table|2| The model spectrum is plotted as the dashed line in Figure 
® 

The maximum brightness temperature of the jet emission does 
not depend on / for model Bl and from equation 1331 we find 
7b,max = 1-1 x 10 ll) K. We can also calculate at what distance / 



along the jet the jet material is dense enough so that relativis- 
tic induced Compton scattering becomes important. From equa- 
tion I34i we find that this loss process would only play a role for 
I < 1 .5 x 10~ 7 which is well inside / m ; n . We therefore conclude that 
our model can be applied to the jet of Cygnus X-l. 

The jet radius ro is very small. The ballistic jet model is coni- 
cal and therefore we can define a jet opening angle as 9 = 2rQ/xo = 
5" which is much smaller than the observational upper limit of 2° 
IStirling et alll200ll) . It is of course possible to assume a larger ra- 
dius for the jet by dropping the assumption of equipartition. How- 
ever, setting the jet radius equal to the observational upper limit 
would imply that the jet material is out of equipartition by several 
orders of magnitude. 

The physical limits of the jet flow, / m ; n and Z max , can be de- 
rived from equation <12> by setting % (/min) = T (/max) = 1 for those 
frequencies at which the flat spectrum is required to break to the op- 
tically thick or thin regime, respectively. In our example of Cygnus 
X-l we used a lower break of just under 100 MHz and an up- 
per break of 3 x 10 14 Hz corresponding to a wavelength of 1 pm. 
The resulting limit Z m ; n = 4.3 x 10~ 5 , associated with the break to 
optical thin conditions, corresponds to a distance of about 6000 
Schwarzschild radii from the central black hole with a mass of 
10 M tGierli hski et allll999h . Since for the ballistic jet model B 1 
/ m j n decreases linearly with the break frequency, the lower physical 
limit decreases to about 500 Schwarzschild radii if the flat spec- 
trum extends to the near-UV. Clearly the determination of the high 
frequency break of the flat part of the spectrum can put interesting 
constraints on the distance from the central black holes at which 
jets become ballistic. Similarly the low frequency break constrains 
the overall extent of the ballistic jet flow. 

The observable extent of the jet flow depends linearly on 
the observing frequency. Since / max (8.4GHz) = 1, we would ex- 
pect that / max (15GHz) ~ 0.6. This prediction significantly exceeds 
the tentative ex tension of the Cygnus X-l jet of / max (15 GHz) < 
0.13 reported in Stirli ng et aljil99 8). However, the observations at 
15 GHz were not simultaneous with those at 8.4 GHz. 

The strength of the magnetic field, Bo, at xq implies a field 
strength of around 50 T at the distance / m ; n . If / lran is reduced be- 
cause the flat spectrum extends to the near-UV, then the magnetic 
field in the jet has a strength of 500 T about 500 Schwarzschild radii 
away from the black hole. Again this demonstrates that our model 
can provide useful constraints on the conditions at the very base of 
the observed jets despite them not being spatially resolved. 

The strength of the magnetic field, Bq, and the constant in the 
expression for the density of relativistic electrons, Ko, can be used 
to estimate the power of the jet. We find that the energy transport 
rate associated with the magnetic field and the relativistic particles 
alone is 5.2 x 10 25 W. A further 1.1 x 10 25 Wis added by the kinetic 
energy of the electrons. Not surprisingly these numbers are compa- 
rable to the estimates using the lBlandford & Konigll Jl979l) model 
iFender et alll200fj) . If there is a cold proton for every relativistic 
electron in the jet, then its energy transport rate in terms of kinetic 
energy is 1.1 x 10 28 W. This power is about one order of magnitude 
below the time-averaged energy transport rate of the Cygnus X-l 
jet, recently estimated from the observed interaction of the jet with 
the surrounding ISM <Gallo et all 120051) . The flux density of the 
flat jet spectrum does not vary significantly between available ob- 
servations. Unless the jet power varies considerably over timescales 
longer than the timespan since the first available radio observations 
of Cygnus X-l, then our estimate strongly suggests the presence of 
a proton-electron plasma in the jet. 
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Figure 5. The maximum extent of the jet for an observing frequency of 
8.4 GHz in the adiabatic jet model A3 as a function of the magnetic field 
strength Bo- The solid line shows l' m . ix for initial equipartition. The dotted 
line indicates l' max = 1 as required by the observations. The hatched region 
bounded by the long-dashed curve contains the allowed combinations of Bo 
and l' max for the assumption that / m j„ is larger than three Schwarzschild radii 
of a 10 Mq black hole and that the flat spectrum extends to the near-IR. The 
short-dashed line gives l' max for an adiabatic jet out of equipartition with 
/ = 10~ 6 . The model parameters given in Tablel3lare derived for the point 
where the short-dashed and dotted lines cross inside the hatched region. 



4.7.2 Adiabatic jet models 

All three adiabatic jet models are consistent with a partially flat 
emission spectrum. In fact, if we choose the geometrical param- 
eter a\ appropriate for a flat spectrum for each model, then the 
differences between the adiabatic models become small as the dif- 
ferences in the behaviour of the magnetic field as a function of / 
between them is compensated for by the different degrees of con- 
finement of the jet. Therefore and to simplify the discussion below, 
we focus on the adiabatic model A3 with an isotropic magnetic 
field. None of the conclusions change greatly for models Al and 
A2. 

For the adiabatic jet model A3 our assumptions with equation 
<37l lead to 

F v ~ 8.8 x lO-^lxorofl-^o 1/2 v 5 / 2 C x 5/3 mjy, (46) 

where we have set a = 1/5 to allow for a flat section in the spec- 
trum. The calculation of y max now requires /thick » which must be de- 
termined from the implicit equation I26t for a given magnetic field 
strength, Sq. Figure|5|demonstrates that the allowed minimum for 
'max exceeds unity for the assumption of initial equipartition. Since 
observations require / max = 1 at 8.4 GHz, the adiabatic models are 
incompatible with observations unless we drop the requirement of 
equipartition. 

We now introduce a reduction factor / such that the initial en- 
ergy density of the relativistic electrons is a fraction / of the initial 
energy density of the magnetic field. An example of the results for 
'max f° r / = 10 is shown in Figure[5] There are now two possible 
solutions for the strength of the magnetic field. However, we also 
require that the lower size limit of the jet, Z mm , accommodates a 
break of the flat spectrum to the optically thin regime in the near- 
IR. This lower limit cannot lie inside the last stable orbit of the cen- 
tral black hole. Thus we obtain another constraint on the solution 
because > where R$ = 3 x 10 4 m is the Schwarzschild ra- 
dius of a 10 Mq black hole. The additional constraint is also plotted 



Figure 6. Spectrum of the adiabatic jet model A3 without initial equipar- 
tition using the parameters in Table l3l The solid line shows the spectrum. 
The dashed line is the optical depths, T m!lx , at / m j„. The dotted lines indicate 
the frequency for which T max = 100. As expected, the peak in the spectrum 
occurs at T max ~ 1 . 



Model parameter 


Value 


xo 


47 AU 


ro 


1.3 x 10 9 m 


V 


2.5 


Bo 


1.2T 


KO 


7.6 x 10 _8 J 15 m- 3 


D 


2 kpc 


An in 


2.0 x 10~ 8 




200 


Ymax (/min) 


10 6 


l 'j 


0.97 c 




40° 



Table 3. Parameters for the adiabatic jet model A3 without initial equipar- 
tition. The model spectrum is shown in Figurelol 



in Figure|5| Only one of the two possible solutions for / = 1CP 6 
is consistent with this constraint. It is also interesting to note that 
the solution requiring initial equipartition is also inconsistent with 
a physical meaningful value of / mm . 

Using the remaining solution for / = 10~ 6 as an example, 
we can compute all remaining model parameters which are sum- 
marised in Table [3] Figure [5] shows the resulting spectrum. It is 
flat over a wide range of frequencies with the required flux level. 
The dip at around 10 14 Hz is a result of the diminishing optical 
depths of the jet to radiation emitted by electrons with the limiting 
Lorentz factor y max = Ythick ('Ahick) _2a • For higher frequencies 



those parts of the jet close to / m i n which still contain electrons with 
Lorentz factors in excess of Ythick contribute to the emission and 
cause the peak. The position of the peak is located at the frequency 
with optical depth T max ~ 1 at / mm . Clearly the spectrum is not flat 
from radio to near-IR frequencies because of the emission peak. A 
flat spectrum extending over the entire radio to near-IR range could 
be achieved by moving the peak to higher frequencies. However, an 
appropriate adjustment of the model parameters would also tighten 
the constraints on the reduction factor / and thereby on Sq. 

The maximum brightness temperature of the adiabatic model 
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for efficient Compton scattering. The distance at which relativistic 
induced Compton scattering would become important in the adia- 
batic jet is / = 1 .7 x 10~ 16 . Again we are justified to neglect energy 
losses due to Compton scattering. 

The jet radius, ro, for the adiabatic jet is larger than for the 
ballistic jet. Formally, we cannot define a jet opening angle for 
adiabatic jets discussed here, because their shape is not conical. 
However, if a conical shape was assumed, then the opening angle 
inferred from the radius at xq would be 1.3'. 

In our example, the lower limit of the physical extent of the 
jet, / m ; n , corresponds to 4.7 R$. While this is close to the theoretical 
limit and gravitational redshift would affect the spectrum, a further 
decreased value of the reduction factor would increase this limit 
at the expense of requiring an even stronger magnetic field. For 
adiabatic jets the constraints on the nature of the jet flow extend to 
even closer distances from the central black hole than in the ballistic 
case. 

For the adiabatic jet model A3 the observable extent of the 
jet flow is proportional to v~ 3// ( 8fl '. In our example this relation 
implies that /' (15 GHz) ~ 0.3, w hich is still larger than the ob- 
served value o f lStirling etall fi9981. but closer than the prediction 
of the ballistic model. However, as mentioned above, the observa- 
tions are not simultaneous and that may explain the discrepancy in 
both cases. 

The strength of the magnetic field of 1 .2 T is high at xq . How- 
ever, because of the much more collimated geometry of the jet, the 
field strength only increases to 140 T at / m j n . Nevertheless, the en- 
ergy transport rate of the jet due to the magnetic field is very large 
with 10 34 W determined at / m ; n . Due to the small value for the re- 
duction factor /, the contribution of the relativistic electrons to any 
energetic considerations is negligible, even if the kinetic energy of 
possibly associated protons is taken into account. The derived jet 
power exceeds by far all prev ious estimates and is inconsistent with 
the time-averaged jet power iGallo et allEool . unless the jet flow 
is suppressed for long periods on very long timescales. Note also, 
that this jet power corresponds to one hundred times the Eddington 
limiting luminosity of a 10 Mq black hole. The large jet power is 
caused by the significant reduction in the radiative efficiency of the 
synchrotron process well away from equipartition conditions. 



5 CONCLUSIONS 

We construct a model for the synchrotron emission of partially self- 
absorbed jets. The model does not invoke a re-acceleration process 
for the relativistic electrons. All electrons are accelerated only once 
at the lower physical limit of the jet, l^. It is not necessary to 
postulate an unknown re-acceleration mechanism as was suggested 
in previous work l Blandford & Konigl, 1979), because synchrotron 
self-absorption counteracts excessive energy losses. 

Two classes of models are considered. The ballistic jets have a 
conical geometry and their contents do not suffer adiabatic energy 
losses. This situation may arise when jets are initially highly over- 
pressured with respect to their environments. They expand unim- 
peded and random thermal energy is converted into bulk kinetic 
energy, but not dissipated to any external medium. At the end of 
this very rapid expansion the mean free path of the jet material ex- 
ceeds the physical dimensions of the jet itself and follows ballistic 
trajectories. 

The other class of models considered here are adiabatic jets 
confined by either the pressure of the gas surrounding them, surface 
magnetic fields or both. Their geometrical shape is dictated by the 



details of the confinement mechanism which is not the subject of 
this paper. The jet material dissipates energy to the external gas 
during its adiabatic expansion. 

Both classes of models can predict flat emission spectra if en- 
ergy losses of individual electrons are neglected. The ballistic jet 
with a magnetic field perpendicular to the jet axis produces a flat 
spectrum without further assumptions. The adiabatic jet models re- 
quire a specific jet geometry to allow for flat emission spectra. 

Both model classes can predict flat emission spectra, even 
when taking energy losses of the electrons and the magnetic field 
into account. Synchrotron self-absorption prevents radiative energy 
losses below a critical Lorentz factor Ythick- In tne ballistic case 
Ythick is constant along the jet flow and because adiabatic energy 
losses are absent, the energy distribution of the relativistic elec- 
trons remains stationary. The emission prope rties of the jet in this 
case a re essentially identical to the model of Blandford & Konigl 
1 1979). For the adiabatic jet the with a perpendicular magnetic field 
electrons with Lorentz factors at and below Ythick continue to lose 
energy because of the sideways expansion of the jet. However, for 
a geometrical shape given by r k x 1 / 5 the resulting spectrum is 
again flat. Other configurations of the magnetic field can also lead 
to flat emission spectra. The spectral slope is very sensitive to the 
jet shape. For example, slightly changing the jet shape to r« x 1 / 4 
results in a spectrum with F v v 0,38 . 

We show an application of the model to observations of a re- 

I III. 
solved jet in the X-ray binary Cygnus X-l l Stirling et al., 200JD. As 

input we use the flux density of the jet in the flat part of the spec- 
trum, the physical extent of the resolved jet and an assumed break 
of the flat spectrum to optically thin conditions in the near-IR. Both 
model classes can be made consistent with the observational con- 
straints. However, in doing so the adiabatic jet models require a sig- 
nificant departure from energy equipartition between the magnetic 
field and the relativistic particles. The associated reduced radiative 
efficiency of the jet plasma implies extremely high energy transport 
rates for the jet of around 10 34 W. This jet power exceeds the Ed- 
dington limiting luminosity of a 10 Mq black hole by two orders of 
magnitude. 

The ballistic jet is consistent with current observations and 
require s energy transport rates well below the time-averaged jet 
power l Gal lo et alll2005h . This result holds even if the kinetic en- 
ergy of one non-radiating proton per relativistic electron is taken 
into account. Unless Cygnus X-l ceases to produce a jet for long 
periods of time, then the ballistic model requires a proton-electron 
jet plasma to explain the large accumulated energy in the region 
where the jet interacts with the surrounding gas. 

The jet opening angle of 5" required by the ballistic model is 
very small. The velocity perpendicular to the jet axis established 
in the jet acceleration region must be very small for this model to 
work. If the ballistic jet results from a rapid initial expansion of a 
highly overpressured jet, then the requirements are even more se- 
vere. This result places stringent limits on the acceleration process. 

Both model classes can probe the conditions in jets on scales 
unresolvable with current telescopes. The most important measure- 
ment in this respect is the high frequency cut-off of the flat spec- 
trum which is formed closest to the jet acceleration region. In or- 
der to resolve the remaining uncertainties with the model, the most 
helpful measurements would be to resolve the jet along its axis in 
quasi-simultaneous observations at two different frequencies. As 
mentioned above, even a factor 2 between the observing frequen- 
cies employed would help to constrain the geometrical shape of the 
jet and thereby help to decide which of the model classes is com- 
patible with observations. 
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Obviously from the work presented here we cannot rule out 
the possibility of continued re-acceleration of the radiating elec- 
trons in the jet. However, such an energetisation process is not 
necessary to produce flat spectra from self-absorbed, synchrotron 
emitting jets. 
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